arXiv:1508.07728vl [hep-th] 31 Aug 2015 


Pis’ma V ZhETF 


INR-TH/2015-022 


Towards conformal cosmology 

M.Libanov+, V.Rubakov+*, 2 ) G.Rubtsov+ 3) 

+ Institute for Nuclear Research of the Russian Academy of Sciences, 

60th October Anniversary Prospect, 7a, 117312 Moscow, Russia 

* Department of Particle Physics and Cosmology, Physics Faculty, Moscow State University 
Vorobjevy Gory, 119991, Moscow, Russia 


Approximate de Sitter symmetry of inflating Universe is responsible for the approximate flatness of the 
power spectrum of scalar perturbations. However, this is not the only option. Another symmetry which can 
explain nearly scale-invariant power spectrum is conformal invariance. We give a short review of models based 
on conformal symmetry which lead to the scale-invariant spectrum of the scalar perturbations. We discuss 
also potentially observable features of these models. 


Observational data show that primordial scalar per¬ 
turbations in the Universe must have been generated at 
some early cosmological stage, preceding the hot epoch. 
They are nearly Gaussian and have nearly flat power 
spectrum [T]. The first property suggests that these 
perturbations originate from amplified vacuum fluctua¬ 
tions of weakly coupled quantum field(s). Indeed, the 
defining property of Gaussian random field C(x) is that 
it obeys the Isserlis-Wick theorem, which holds also for 
any free quantum field in its vacuum state, while lin¬ 
ear evolution in classical background does not induce 
non-Gaussianity. 

The second property is also very suggestive. The 
power spectrum V(k) defined as 

(C(k)C(k') = ^P(fc)<5(k + k') 
gives the fluctuation in logarithmic interval of momenta, 
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where Ug is a spectral index. The flat or scale invariant 
spectrum corresponds to Ug = 1. In early 70’s Harrison, 
Zeldovich and Peebles and Yu [5] conjectured that the 
spectrum is flat, to avoid large deviation from homo¬ 
geneity and isotropy of the observed Universe on large 
scales and black hole formation on small scales. Gurrent 
observational data [T] show that the power spectrum is 
indded nearly flat and give rig — I ~ —0.032. 

The flatness of the power spectrum may be due 
to some symmetry. The best known candidate is 
the symmetry SO{A, I) of the de Sitter metric ds^ = 
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dt^ — e^^*c?x^, which includes spatial dilatations sup¬ 
plemented by time translations: x -A Ax, t ^ t — 
(H)~^\ogX. This is the approximate symmetry of in¬ 
flating Universe [3], and, indeed, the inflationary mech¬ 
anism of the generation of scalar perturbations [3] pro¬ 
duces almost flat power spectrum. Despite this success 
of the inflationary paradigm, search for alternatives to 
inflation in general and to de Sitter symmetry in par¬ 
ticular is of obvious interest. Alternatives to inflation 
include ekpyrotic model with negative exponential 
potential and bounce [6], ’’starting” [7] Universe, etc, 
while there are several mechanisms capable of producing 
flat or almost flat scalar spectrum [iinmniin]. In some 
cases, there is no obvious symmetry that guarantees the 
flatness, i.e., the scalar spectrum is flat accidentally. 

In quest for an alternative symmetry behind the 
nearly flat scalar spectrum one naturally turns to con¬ 
formal symmetry 5'0(4, 2) [T^ |T31 HH US] (see Ref. m 
for related discussion). Gonformal group includes di¬ 
latations, -A Xx^, which in the end may be responsi¬ 
ble for the scale-invariant scalar spectrum. An assump¬ 
tion of conformal invariance at the time the primordial 
perturbations are generated is in line with the viewpoint 
that the underlying theory of Nature may have confor¬ 
mal phase, and that the Universe may have started off 
from, or passed through an unstable conformal state 
and then evolved to much less symmetric state we see 
today. 

In this paper we give a short review of models 
based on conformal symmetry, which lead to the scale- 
invariant spectrum of the scalar perturbations. It is 
worth noting that one does not necessarily have to con¬ 
sider these models as alternatives to inflation, since 
some of them can work in inflating Universe as well. 
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1. TWO SAMPLE WAYS OF GETTING FLAT 
SCALAR SPECTRUM 

11 Conformal and global symmetry instead of 
de Sitter symmetry 

To begin with, let us consider a scenario proposed in 
Ref. |12) . In this scenario conformal symmetry is sup¬ 
plemented by a global symmetry. The simplest model 
of this sort has global symmetry 17(1) and involves com¬ 
plex scalar field cj), which is conformally coupled to grav¬ 
ity. The action of the model is S' = Sq+m + S^, where 
Sc+M is the action of gravity and dominating matter, 
while the non-trivial dynamics of the scalar field, which 
is assumed to be spectator, is governed by 

Thus, quartic potential allowed by conformal invariance 
is assumed to be negative. Therefore, (() = 0 is an unsta¬ 
ble state with unbroken conformal symmetry. One as¬ 
sumes that the background space-time is homogeneous, 
isotropic and spatially flat, ds^ = a?{rii){drf — dx^). 
Then in terms of the field y(77,x) = a{r])(j)(r],x) the 
dynamics is the same as in flat space-time. 


An important assumption of the entire scenario is that 
the rolling stage begins early enough, so that there is 
time at which the following inequality holds: 

k{T]^ - ?7) » 1, (4) 

where k = |k| is conformal momentum. Since the mo¬ 
menta k of cosmological significance are as small as the 
present Hubble parameter, this inequality means that 
the duration of the rolling stage in conformal time is 
longer than the conformal time elapsed from, say, the 
beginning of the hot Big Bang expansion to the present 
epoch. This is only possible if the hot Big Bang stage 
was preceded by some other epoch, at which the stan¬ 
dard horizon problem is solved; the mechanism we dis¬ 
cuss here is meant to operate at that epoch. We note 
in passing that the latter property is inherent in most, 
if not all, mechanisms of the generation of cosmological 
perturbations. 

Equation (|31) is exactly the same as equation for min¬ 
imally coupled massless scalar field in the de Sitter back¬ 
ground. Nevertheless, let us briefly discuss its solutions. 
At early times, when the inequality (|1]) is satisfied, the 
third term in © is negligible and 6 x 2 is free massless 
quantum field. 


84 ,= I d^Xy/^ 


= ( 1 ) 

Spatially homogeneous background approaches the late¬ 
time attractor 


SX 2 {^,V)=J ^ 2 ^^y 2 ^ ^(k,x,ry)ik + h.c.) , 
whose modes are 


Xciv) 


1 

Hv* - v) ’ 


( 2 ) 


where 77 * is an arbitrary real parameter, and we con¬ 
sider real solution, without loss of generality. The back¬ 
ground solution m breaks conformal group 5'0(4, 2) —>■ 
SO{A, 1). The meaning of the parameter 77 * is that the 
field Xc would run away to infinity as 77 —>■ 77 *, if the 
scalar potential remained negative quartic at arbitrar¬ 
ily large fields. It is worth noting that the particular 
behaviour Xc oc ( 77 * — 77 )“^ is dictated by conformal sym¬ 
metry. 


Phase perturbations. To see how the scale invari¬ 
ant spectrum emerges in the model, let us consider 
perturbations of the phase 9 = Aigcj), or, for the real 
background perturbations of the imaginary part 
X 2 = Imx/V^. At the linearized level, perturbations 
of the phase and modulus decouple and the linearized 
equation is 


dx^^(k,x,77)=e*“V (5) 

Here Ak and are annihilation and creation op¬ 
erators obeying the standard commutational relation, 
[Ak, a]).,] = d(k — k'). It is natural to assume that the 
field ( 5 x 2 is initially in its vacuum state. 

The rolling background Xciv) produces an effective 
“horizon” for the perturbations 6 x 2 - The oscillations 
© terminate when the mode exits the “horizon”, i.e., 
at k{r]^ — rf) ~ 1. The solution to eq. (jS]) with the initial 
condition ([SI) is 

(5x^^ (k, X, 77 ) = e*kx-.7cr,. . ( 6 ) 

where 

= (7) 

and is the Hankel function. In the late-time super- 
”horizon” regime, when 7 ( 77 * — 77 ) <C 1 and the third 
term in © dominates, one has 


<5X2 = 0 . (3) 


F(fc, 77 * - 77 ) = 


(<5X2)" - d,d, 6x2 


( 77 * - 77)2 


Hv* - 9) ' 


( 8 ) 
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Hence, the super-” horizon” perturbations of the phase 
S0 = 5 x 2 1Xc are time-independent, 


(56*(x) 


<5X2 (x, 7?) 


d^k 

47 ^ 3 / 2 ^ 3/2 




(9) 

This expression describes Gaussian random field whose 
power spectrum is flat: 




{ 2 ny 


( 10 ) 


We emphasize that this result is an automatic con¬ 
sequence of the global 17(1) and conformal symmetries. 
To see this, let us consider long wavelength regime. In 
this regime the second term in ([3]) is negligible and ([3]) 
becomes the equation for spatially homogeneous pertur¬ 
bation. Recall that Xc is the spatially homogeneous so¬ 
lution to the full field equation dD, hence, due to 17(1) 
symmetry, e*“xc is also a solution, where a is a real 
constant. For small a the latter solution is Xc + *Q<Xc) 
and the imaginary part is a small perturbation, which 
is precisely 5 x 2 - So, if the perturbation 6 x 2 oscillates 
with unit amplitude, it behaves at late times as (cf. ([5]) 
and dH])) 5 x2 = C/[k{r]^, — rj)], where the factor k~^ is 
evident on dimensional grounds and C is independent 
of time and k. 

Deviation from exact conformal invariance naturally 
gives rise to the tilt in the power spectrum, which de¬ 
pends both on the way conformal invariance is broken 
and on the evolution of the scale factor m- 


12 Galilean Genesis 

The Galileon model has been introduced in Ref. [18] . 
In Minkowski space-time, the Lagrangian of the simplest 
conformally-invariant version m of the model is 


= f‘^e^^d^TTd>"TT + j^df.TTd'^TTnn + , 

( 11 ) 

where □ = . This Lagrangian is conformally invari¬ 

ant, with TT transforming under dilatations as 


In the Galilean Genesis scenario the Universe be¬ 
gins from Minkowski space-time. The field equation in 
Minkowski space-time admits a homogeneous attractor 
solution 


Hcit* — t) 


( 12 ) 


where Hq = 2A^/(3/). The form of the solution is 
again dictated by conformal invariance. 


Initial energy density is zero, while effective pressure 
is negative. Then the energy density slowly builds up 
and the Hubble parameter grows in time, 

until (G — t)^ Hq^ ■ f /Mpi. The growth of the Hubble 
parameter is due to violation of all energy conditions. 
Nevertheless, the theory is fully self-consistent: there 
are no ghosts, tachyons, and other pathologies (there 
is, however, superluminality issue which has not been 
quite settled [19]). At some time Galileon is assumed 
to transmit its energy to conventional matter, and hot 
epoch begins. 

Galileon perturbations per se are not suitable for 
generating scalar perturbations (see discussion in Sec.Q. 
For the purpose of generating the scalar perturbations 
another field 0 of conformal weight 0 is introduced. By 
conformal invariance, its quadratic Lagrangian has the 
form 

Le=e^-{d,eY ^ = ( 9 ^ 0 ) 2 . 

That is, the dynamics of perturbations 59 in the back¬ 
ground TTc is exactly the same as in the conformal rolling 
model discussed above. 

The similarity between Galilean Genesis and confor¬ 
mal rolling model is not accidental. In Ref. [3D] general 
arguments are given, which show that the scale invariant 
power spectrum is inherent in an entire class of models. 
The general setting is conformally invariant theory of a 
scalar field p of conformal weight A 7 ^ 0 in effectively 
Minkowski space-time. Up to rescaling this field corre¬ 
sponds to |(/>| in the conformal rolling model and to e'"’ 
in the Galilean Genesis scenario; in both models A = 1. 
The form of homogeneous classical solution 

is dictated by conformal invariance (t is conformal time 
in the case of the conformal rolling scenario). As men¬ 
tioned above, the perturbations of the field p do not have 
flat power spectrum, so another spectator scalar field 9 
of conformal weight 0 is introduced. Then by confor¬ 
mal invariance, the kinetic term in the Lagrangian of 
9 is Lg (X p^/^{d^9Y. Assuming that possible poten¬ 
tial terms are negligible the Lagrangian in the rolling 
background (1131) takes the form 


Le 


const 
{t* - ty 


{d,9f, 


which is exactly the same as the Lagrangian of a scalar 
field minimally coupled to gravity in de Sitter space with 
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conformal time t and scale factor a{t) oc !/(<, — t). As 
a result, 9 develops perturbations with the flat power 
spectrum. 

13 Further aspects 

Obviously, generating the field perturbations S9 is 
not the whole story. There are several other ingredients 
of the conformal scenario. Some of them have not yet 
been worked out in detail. 

Beginning of rolling. The rolling stage (USD has to 
start in one or another way. One possibility is sponta¬ 
neous decay of the unstable conformally invariant vac¬ 
uum p = 0, which proceeds through bubble nucleation. 
Such a decay has been discussed in Ref. [21] within the 
holographic approach (holographic picture of the confor¬ 
mal scenario has been suggested earlier in Refs, 

Even though the rolling field is not spatially homoge¬ 
neous in the false vacuum decay process, the field per¬ 
turbations have the properties we discuss in Secs. HT] 
0 . 



Fig. 1. The scalar potential in the pseudo-Nambu- 
Goldstone scenario. 

erty of the conformal scenario, however, as it holds in 
all models employing the curvaton or modulated decay 
mechanism, including versions of inflation. 


End of rolling. The rolling stage (USD should termi¬ 
nate at some late time. This implies that conformal 
invariance is broken at large field values. As an exam¬ 
ple, in the conformal rolling scenario of Sec. [TTjone as¬ 
sumes that the scalar potential has a minimum or nearly 
flat valley at large \(j)\ and that \(j)\ eventually settles 
there. Furthermore, in the Galilean Genesis scenario 
there must be a stage of defrosting, i.e., transmission of 
energy from the rolling field p to heat [23], after which 
the usual hot Big Bang epoch begins. 

Reprocessing perturbations 59 into adiabatic 
perturbations. The field perturbations 59 are to be 
converted into adiabatic perturbations. This can hap¬ 
pen at the hot Big Bang epoch. One possibility is to 
make use of the curvaton mechanism ESI- As an exam¬ 
ple, in the scenario of Sec. [TT] the phase 9 may actually 
be a pseudo-Nambu-Goldstone held. Generically, con¬ 
formal rolling ends up at a slope of its potential, see 
Fig. [TJ The perturbations 59 are reprocessed into adi¬ 
abatic perturbations at later time at the hot Big Bang 
epoch when the held 9 oscillates and decays into con¬ 
ventional particles, cf. Ref. ^5] . 

Another possibility is the modulated decay mecha¬ 
nism suggested in the inhationary context in Refs. [221 

EHmSj. 

In both cases the conversion of the held perturba¬ 
tions 59 into adiabatic perturbations induces some de¬ 
gree of non-Gaussianity. This is not a specihc prop- 


2. PERTURBATIONS OF MODULUS 


Let us now come back to the conformal rolling stage. 
From now on we use the nomenclature of the negative 
quartic model of Sec. [11] for dehniteness. We are inter¬ 
ested in the radial perturbations [Hisniisi] , i.e., per¬ 
turbations of the modulus of the held x, or, with our 
convention of real background Xc perturbations of the 
real part xi = Re x/v^- At the linearized level, they 
obey the following equation, 

{5xi)" - d^di 5xi - 6/i^Xc^Xi = 0 ■ 


Its solution that tends to properly normalized mode of 
free quantum held at early times, k{r]^: —??)—>■ oo, is 


,5x1 = — 

Att 




77)]-Bk+h.c. , 


where i?k, is another set of annihilation and cre¬ 
ation operators. At late times, when fc(? 7 * — 77 ) <C 1 
(super-”horizon” regime), one has 


5xi = 


3 1 

471^/2 




■ Rk + h.c. . 


Hence, the super-”horizon” perturbations of the modu¬ 
lus have red power spectrum 


'P\^\{k) (xk (14) 

The dependence 5xi oc (ry* — is interpreted in 
terms of the local shift of the “end time” parameter 77 *. 
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Indeed, with the background field given by ©, the sum 
Xc + <Jxi, i-e., the radial field including perturbations, is 
the linearized form of 

where 77 * (x) = 77 * + 6t]^ (x) and 


%(*) = -4^ / ^ ■ 

(16) 

So, the infrared radial modes modify the effective back¬ 
ground by transforming the “end time” parameter 77 * 
into time-independent random field that slowly varies in 
space, with red power spectrum clearly seen from dT51) . 
This observation is valid beyond the linear approxima¬ 
tion: once the spatial scale of variation of xi(x, 77 ) ex¬ 
ceeds the “horizon” size, spatial gradients in eq. (HD are 
negligible, and the late-time solutions to the full non¬ 
linear field equation have locally one and the same form 
(ED, modulo slow variation of 77 * in space. 

A few remarks are in order. First, the infrared modes 
contribute both to the field ^ 77 * (x) itself and to its spa¬ 
tial derivatives. The contribution of the modes which 
are superhorizon today, i.e., have momenta k < Hq, to 
the fluctuation of diC]^ is given by 


(a,774x)5j774x))fc<^ f ^ 

47 r Jk<Ho * 


. 3h2 Ho 

= ^*A^log— , (17) 


where A is the infrared cutoff which parametrizes our 
ignorance of the dynamics at the beginning of the con¬ 
formal rolling stage. 

Second, modulo field redefinition and notations, the 
properties of Galileon perturbations are exactly the 
same as the properties of radial perturbations 5xi in 
conformal rolling scenario [32]. Furthermore, these 
properties are unambiguously determined by conformal 
invariance |3H [50] . The same properties - fiat and red 
spectra of zero conformal weight and rolling fields, re¬ 
spectively - are inherent in the false vacuum decay setup 
mentioned in Sec. Thus, we are dealing with the 
whole class of models. 

Finally, so far we have considered spectator fields. 
That is, we have ignored the backreaction of the scalar 
fields on gravity. In particular, we have considered a 
scalar field conformally coupled to gravity (in the con¬ 
formal rolling scenario of Sec. ED and assumed negligi¬ 
ble energy density (in both Galilean Genesis and con¬ 
formal rolling cases). However, it is of interest to con¬ 
sider also dynamical versions with scalar fields mini- 



Fig. 2. Due to the perturbations of the radial field, the 
evolution of phase perturbations proceeds in inhomoge¬ 
neous background that slowly varies in space. 


mally coupled to gravity and dominating the cosmolog¬ 
ical evolution. In that case there is a potential danger 
that the strong-coupling regime arises at too low energy 
scales M- This option has been studied in Ref. [33] . 
It has been shown that mixing of the scalar field(s) 
with the metric in dynamical pseudo-conformal models 
does not introduce new strong-coupling UV scales. Fur¬ 
thermore, the spectator approximation gives correct re¬ 
sults in dynamical models provided that the background 
space-time is sufficiently fiat . This applies, in particu¬ 
lar, to potentially observable effects discussed in Sec. 3. 
These effects are inherent in the entire class of both 
spectator and dynamical (pseudo-)conformal models. 


3. EFFECT OF INFRARED RADIAL MODES 
ON PERTURBATIONS OF PHASE 

Let us discuss how the interaction with the infrared 
radial modes affects the properties of the phase pertur¬ 
bations 60 [3nil3T||34]. To this end, we consider pertur¬ 
bations of the imaginary part < 5 x 2 , whose wavelengths 
are much smaller than the scale of the spatial variation 
of the modulus. Because of the separation of scales, 
perturbations 6 x 2 can still be treated in the linear ap¬ 
proximation, but now in the background (fT51) . see Fig.j^j 

Since our concern is the infrared part of 77 , (x), we 
make use of the spatial gradient expansion, consider a 
region near the origin and write 

77*(x) = 77*(0)-7;ia;j-|-... , Ui =-i9j77»(x)|x=o , (18) 

where dots denote higher order terms in x. Impor¬ 
tantly, the field didjCj^iyc) has blue power spectrum, un¬ 
like 77 * (x) and 9 ^ 77 * (x), so the major effect of the infrared 
modes is accounted for by considering the two terms of 
the gradient expansion written explicitly in (1181) . Fur- 
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thermore, we assume in what follows that |v| <C 1. The 
expansion in |v| is legitimate, since the field v(x) has 
fiat power spectrum (cf. (IT71) '). so the fiuctuation of v 
is of order h'^\ log A|, where A is the infrared cutoff, and 
it is small for small h and not too large | log A|. 

Keeping the two terms in only, we have, instead 

of ([HI), 

{ 6 x 2 )” - dA 6 x 2 - r ^ - ^ 5 x 2 = 0 . (19) 

[ 77 * ( 0 ) - vx - 77 ]^ 

We observe that the denominator in the expression for 
the background field 

Ac = , r ^ -T (20) 

/i[77*(0) - 77 - vx] 

contains the combination 77 * (0) — (77 + vx). We interpret 
this as the local time shift and Lorentz boost of the origi¬ 
nal background Note that the field (EHl) is a solution 
to the field equation m in our approximation. Our in¬ 
terpretation makes it clear that the solutions to eq. m 
can be obtained by time translation and Lorentz boost 
of the original solution ([B]), ([7]). In particular, instead 
of ([B]) the phase field freezes out at 


(50(x) = ih 


d^k 


^zkx 


47r3/2\/fc(fc -1- kv) 


(x) ik +/i-c. 


( 21 ) 

So far we discussed the dynamics of the phase per¬ 
turbations 60 at the conformal rolling stage which is 
governed solely be their interaction with the background 
field ( 1201 ) as well as with the radial perturbations Sxi; 
the evolution of the scale factor 0 ( 77 ) is irrelevant. Af¬ 
ter the end of conformal rolling, the situation is re¬ 
versed. Once the radial field jc^j has relaxed to the min¬ 
imum of the scalar potential, the phase 0 is a massless 
scalar field minimally coupled to gravity (this is true 
for any Nambu-Goldstone field [35]). Since we are talk¬ 
ing about a yet unknown pre-hot epoch, it is legitimate 
to ask what happens to the perturbations of the phase 
right after the end of conformal rolling. Barring fine 
tuning, there are two possibilities for the perturbations 
60: 


(i) they are already superhorizon in the conventional 
sense at that time, or 

(ii) they are still subhorizon. 


31 Superhorizon phase perturbations 

Let us consider the first sub-scenario 1301131]: phase 
perturbations do not evolve after the end of the con¬ 
formal rolling stage, and the properties of the adiabatic 
perturbations are determined entirely by the dynamics 


at conformal rolling (modulo possible non-Gaussianity 
generated at the conversion epoch, see Sec. mi). This 
option is particularly natural in the Galilean Genesis, 
but it is not contrived in the conformal rolling scenario 
of Sec. [TT] either. 

In that case, there are no potentially observable ef¬ 
fects to the linear order in v. Indeed, using m one 
finds for, e.g., two-point correlation function 


(60(x)60(x')} oc 


d^k 1 
k (^k + kv)"^ 
d^q 1 
q q^ 


ik{x — 


— (x) — 77 * {x')) 


( 22 ) 


where we have used m in the exponent and changed 
the integration variable from k to 


q = k -I- A:v , g = jqj = fc -|- kv, (23) 


which is nothing but Lorentz boost. The result (EU) is 
precisely the two-point correlation function of the lin¬ 
ear field dU). The latter argument is straightforwardly 
generalized to multiple correlators: for a given realiza¬ 
tion of the random field 77 * (x), they are all expressed in 
terms of the two-point correlation function (1^^ . The 
reason for the disappearance of the linear order effect is 
obviously Lorentz invariance. 

The non-trivial effect of the large wavelength pertur¬ 
bations 6 r]^{x) on the perturbations of the phase, and 
hence on the resulting adiabatic perturbations, occurs 
for the first time at the second order in the gradient 
expansion, i.e., at the order didjVj^: [30]. Let us concen¬ 
trate for the moment on the effect of the modes of 6 'q^, 
whose present wavelengths exceed the present Hubble 
size. We are dealing with one realization of the random 
field < 577 *, hence at the second order of the gradient ex¬ 
pansion, didjr]^ is merely a tensor, constant throughout 
the visible Universe. In this long wavelength regime the 
perturbation of the phase has the following form 


60{k) 


Am^l'^'Jkq 



TT kikj 



* T h.c.j 


(24) 

where q is given by (1231) . It results in the power spec¬ 
trum of the adiabatic perturbation which depends on 
directionality of momentum (see [S^ for details) 


/ ^ 

Vcik) = 'P°{k)il + ci-h--j^ ■ {nk)i{nk)jWij 


. 7,2 . ,.\2 


In the first non-trivial term, Wij is a traceless symmet¬ 
ric tensor of a general form with unit normalization, 
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WijWij = 1, nk is a unit vector, nk = k/A:, and ci is a 
constant of order 1 whose actual value is undetermined 
because of the cosmic variance. The last term is the 
result of the expansion of 69 in h; the corresponding 
term is not present in (12411 . In this term, u is a unit 
vector independent of rcy , and the positive parameter 
C 2 is logarithmically enhanced due to the infrared ef¬ 
fects, C 2 oc log-^. This is the first place where the 
deep infrared modes show up. Clearly, their effect is 
subdominant for small h. 

We see that the large wavelength modes induce sta¬ 
tistical anisotropy in the adiabatic perturbations. The 
statistical anisotropy encoded in the second term in 
(113 is similar to that commonly discussed in inflation¬ 
ary context [30] (see Ref. [37] for earlier analysis), and, 
indeed, generated in some concrete inflationary mod¬ 
els m- it does not decay as momentum increases and 
has special tensorial form (nku)^ with constant u. On 
the other hand, the first non-trivial term in (1251) has the 
general tensorial structure and decreases with momen¬ 
tum. The latter property is somewhat similar to the 
situation that occurs in cosmological models with the 
anisotropic expansion before inflation [39] . 

Surprisingly, the statistical anisotropy in the form of 
the special type quadrupole was found in the WMAP 
5 and 7 years data gni mi m ia. It was argued, 
however, that the anomaly may result from the detec¬ 
tor beam asymmetry not accounted for in the WMAP 
analysis [H]. In the final 9 years data release, WMAP 
collaboration provided a set of maps deconvolved with 
the instrument response function corresponding to the 
beam asymmetry effect gS]. Deconvolved maps do 
not indicate the statistical anisotropy allowing to con¬ 
strain the coupling constant in the first sub-scenario: 
Ai^log-^ < 1.2 g0]. Later, the Planck data led to 
stronger constraints on the statistical anisotropy mm 
and correspondingly to tighter limits on the coupling 
constant, h'^ log < 0.30 [33] . 

Another effect that emerges at order h? is the non- 
Gaussianity of the perturbations 59^ and hence the adi¬ 
abatic perturbations over and beyond the non- 

Gaussianity that may be generated at the time when 
the phase perturbations get reprocessed into the adia¬ 
batic perturbations. In the absence of the cubic self¬ 
interaction of the field 0, the intrinsic bispectrum van¬ 
ishes, so we have to consider the trispectrum. It is fully 
calculated [33], the most striking feature being the sin¬ 
gularity in the limit where two momenta are equal in ab¬ 
solute value and have opposite directions (folded limit, 
in nomenclature of Ref. my- 


(CkiCk 2 Ck 3 Ck 4 ) = const • (5 I 

\i=l 



ki2 — ki k2 —^ 0 , 


i.e., the trispectrum blows up as ki 2 - This is in con¬ 
trast to trispectra obtained in single-field inflationary 
models; indeed, there are general arguments m show¬ 
ing that in these models, the four-point function is finite 
in the limit ki2 —)■ 0. The singularity in the four-point 
function (|26|1 is due to the enhancement of the radial 
perturbations 5\(l)\ at low momenta. 

Even though the results (1331) , (1301) were first derived 
in the concrete model with negative quartic potential 
(Sec.[TT]), they are actually consequences of consistency 
relations m valid in the whole class of conformal mod¬ 
els. Furthermore, these consistency relations enable one 
to calculate the one-loop non-Gaussianity in the folded 
limit ki 2 —>■ 0. Interestingly the one-loop contribution 
to the trispectrum is even more singular in the folded 
limit than the tree-level result ([301): one finds [52] 

(Ck4Ck2Ck3Ck4)^^-'°°^) =constA 

with suppressed coefficient (the suppression factor is 
in the model of Sec. [TT]). The fcC 2 ^-enhancement of the 
one-loop trispectrum in the folded limit makes the non- 
Gaussianity even more promising from the observational 
viewpoint. 

32 Subhorizon phase perturbations 

Let us consider another option: assume that there is 
a long enough period of time after the end of conformal 
rolling, at which the phase perturbations remain sub¬ 
horizon in the conventional sense |34j . This option is 
fairly natural in the conformal rolling model of Sec. [TT] 
and more contrived in Galilean Genesis. 

The behavior of 69 between the end of conformal 
rolling and horizon exit depends strongly on the evolu¬ 
tion of the scale factor at this intermediate stage. In 
order that the flat power spectrum m be not grossly 
modified at this epoch, the scale factor should evolve 
in such a way that the dynamics of 69 is effectively 
nearly Minkowskian. Although this requirement sounds 
prohibitively restrictive, it is obeyed in the bouncing 
Universe, with matter at the contracting stage having 
super-stiff equation of state, p ^ p. It is worth noting 
in this regard that stiff equation of state is preferred at 


1 

ki2kfk^ 
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Fig. 3. Due to the perturbations of the radial field, the 
evolution of phase perturbations proceeds in inhomoge¬ 
neous background. Perturbations S9 oscillate in time at 
early stage (region I), freeze out at time t] = rjx (x) and 
temporarily stay constant (region II) until the end of 
conformal rolling that occurs at t] — » 7 *(x). Then they 
evolve again, now in nearly Minkowskian regime (region 
III), until the horizon exit time rji. Later on (region IV), 
perturbations S9 are superhorizon and stay constant. 


the contracting stage for other reasons [53l [54] and is 
inherent, e.g., in a scalar field theory with negative ex¬ 
ponential potential, like in the ekpyrotic model It is 
known |55j that in models with super-stiff matter at con¬ 
tracting stage, the resulting power spectrum of scalar 
perturbations is almost the same as that of massless 
scalar field in Minkowski space, P(fc) oc In tractable 
bouncing models like those of Ref. [5], the phase per¬ 
turbations evolve almost like in Minkowski space, exit 
the horizon at the contracting stage, pass through the 
bounce unaffected (cf. Ref. [5]), remain superhorizon 
early at the hot expansion epoch and get reprocessed 
into adiabatic perturbations, as discussed in Sec. |T3| 

Let the effectively Minkowskian stage ends at some 
time rji (see Fig. Ej). The field (50(x, 77 *), determined by 
the dynamics at the conformal rolling stage, serves as 
the initial condition for further Minkowskian evolution 
from 77 * to r]i. We are interested in the properties of the 
phase perturbations at ry = 771 , as these properties are 
inherited by the adiabatic perturbations. 

To the leading order in h, we find nothing new: the 
phase perturbations at 77 = 771 are Gaussian and have 
flat power spectrum. Subleading orders in h are more 
interesting. The effect of the perturbations < 577 * (x) on 
the phase perturbations SO is twofold. First, the pertur¬ 
bations Srjif modify the dynamics of SO at the conformal 


rolling stage. This property is the same as in Sec. |5T| 
Second, the initial condition for the Minkwskian evo¬ 
lution is now imposed at the non-trivial hypersurface 
77 = 77 *(x). This is illustrated in Fig. E) The net re¬ 
sult is that the perturbation S9(x) at the time 771 is a 
combination of two Gaussian random fields originating 
from vacuum fluctuations of the phase 0 and radial field 
\(j)\, respectively. This leads to potentially observable 
effects. 

Let us consider the phase perturbation of given mo¬ 
mentum k. At the end of conformal rolling it is given by 
(EU with 856/drj = 0. This gives the initial condition 
for the evolution at the intermediate stage, when the 
perturbation is a linear combination of exp( 7 kx ± ikrj). 
So, after the second freeze-out it is a linear combination 
of waves coming from the direction nk = k/fc and from 
the opposite direction and traveling distance r = 771 — 77 *. 
This leaves an imprint on (50(k) of the random field v 
existing at points x = intT”. In particular, one finds 
the power spectrum with the non-trivial dependence on 
nki 


Vse{ki) = T’o (1 + Hk • [v(x = -l-nkr) - v(x = -nkr)j) . 


As a result, the statistical anisotropy of the adiabatic 
perturbations has the form 

Vc{k)=vf{k)[l + Q{ja^)], (27) 

where is independent of the directionality of mo¬ 
mentum (nearly flat spectrum with small tilt) and 
( 5 (nk) is itself a random field, which depends on the 
direction of k only. Unlike the statistical anisotropy 
discussed in the inflationary context (36] ES |39] and 
also in Sec. ED the function ( 5 (nk) contains all even 
angular harmonics, starting from quadrupole. We give 
here the expression for ( 5 (nk) which accounts for the 
quadrupole component only (see Ref. |34] for all multi¬ 
poles), ( 5 (nk) = where Wij is a general sym¬ 

metric traceless tensor normalized to unity, WijWij = 1 , 
and the variance of the quadrupole component (in the 
sense of an ensemble of universes) is 


(Q^) 


225 / 1 ^ 
327r2 ■ 


(28) 


Of course, the precise values of the multipoles of (3(nk) 
in our patch of the Universe are undetermined because 
of the cosmic variance. It is worth noting also that all 
multipoles are independent of k and, hence, unlike in the 
version of Sec. ED there is no suppression of the leading 
order effect on cosmic microwave background (CMB) 
power spectrum at large 1. This property enables one 
to utilize the high statistics of the Planck data up to 
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I = 1600 and place a strong limit on the coupling con¬ 
stant in the sub-scenario with intermediate stage [49]: 

< 0.0011 . (29) 


The statistical anisotropy is probably the most promis¬ 
ing signature of this sub-scenario. 

The second effect is non-Gaussianity. While, as be¬ 
fore, the bispectrum vanishes, the four-point correlation 
function has a peculiar form 


(CkCkCk'Ck') 


V^°\k) vf\k') 

47rfc3 47rA;'^ 


(5(k-y k)5(k'-P k') 


X [1 + ^^(^(nk • nk')] + (k -O' k^) -|- (k -o- k') . 

(30) 


The leading term in (|30ll (unity in square brackets) is 
the Gaussian part, while the non-Gaussianity is encoded 
in Fmg = 0{h^). Note that the structure of the non- 
Gaussian part is fairly similar to that of the discon¬ 
nected four-point function. Note also that F^g depends 
on the angle between k and k' only. If the angle between 
k' and k is small, i.e., |nk — nk'| <C 1 , the leading be¬ 
havior of Fmg is 


Fng 


“ih? 

^log 

TT^ 


const 

nk - nk' 


where constant in the argument of logarithm cannot be 
reliably calculated because of the cosmic variance. The 
logarithmic behavior does not hold for arbitrarily small 
I nk — nk' I: the function (nk — nk') flattens out most 

likely at |nk — nk'| [k{r]i — and certainly at 

|nk — nk'l [k{r]i — ??*)]“^. So, the parameter (771 — ry*) 
is detectable in principle. 

The third effect is negative scalar tilt 


However, this is not a particularly strong result, as small 
scalar tilt in our scenario may also originate from ex¬ 
plicit violation of conformal invariance at the conformal 
rolling stage [T7] and/or not exactly Minkowskian evo¬ 
lution of 59 at the intermediate stage. Moreover, to 
account for the whole scalar tilt detected by WMAP 
and Planck, one needs h ~ 0.6, in conflict with the con¬ 
straint (ESI). 


4. CONCLUSIONS 

Flat or nearly flat power spectrum of the adia¬ 
batic perturbations may be a consequence of confor¬ 
mal symmetry rather than de Sitter symmetry. Models 


of this sort include conformally coupled complex scalar 
field with negative quartic potential, Galilean Gene¬ 
sis and decay of conformally invariant metastable vac¬ 
uum. Properties of the perturbations in these models 
are to large extent dictated by conformal invariance and 
the predictions are mostly model-independent, at least 
at the leading non-linear level (modulo effects due to 
conversion of field fluctuations into adiabatic perturba¬ 
tions). A peculiar property which has potentially ob¬ 
servable consequences is fluctuations along rolling di¬ 
rection. These fluctuations have red power spectrum 
and can be interpreted in terms of the local time shift. 
Interplay between phase perturbations, responsible for 
density perturbations in the end, and local time shift 
yields the non-trivial correlation properties of the den¬ 
sity perturbations such as statistical anisotropy and the 
intrinsic non-Gaussianity of special forms. The latter 
properties are potentially observable with CMB data. 

This work has been supported by Russian Science 
Foundation Grant No 14-12-01430. 
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